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ABSTRACT 

"Transitional" T Tauri disks have optically thin holes with radii > 10 AU, yet accrete up to the 
median T Tauri rate. Multiple planets inside the hole can torque the gas to high radial speeds over 
large distances, reducing the local surface density while maintaining accretion. Thus multi-planet 
systems, together with reductions in disk opacity due to grain growth, can explain how holes can be 
simultaneously transparent and accreting. There remains the problem of how outer disk gas diffuses 
into the hole. Here it has been proposed that the magnetorotational instability (MRI) erodes disk 
surface layers ionized by stellar X-rays. In contrast to previous work, we find that the extent to which 
surface layers are MRI-active is limited not by ohmic dissipation but by ambipolar diffusion, the latter 
measured by Am: the number of times a neutral hydrogen molecule collides with ions in a dynamical 
time. Simulations by Hawley & Stone showed that Am ~ 100 is necessary for ions to drive MRI 
turbulence in neutral gas. We calculate that in X-ray-irradiated surface layers, Am typically varies 
from ^10~ 3 to 1, depending on the abundance of charge-adsorbing polycyclic aromatic hydrocarbons, 
whose properties we infer from Spitzer observations. We conclude that ionization of H2 by X-rays and 
cosmic rays can sustain, at most, only weak MRI turbulence in surface layers 1-10 g/cm 2 thick, and 
that accretion rates in such layers are too small compared to observed accretion rates for the majority 
of disks. 

Subject headings: accretion, accretion disks — instabilities — ISM: molecules — MHD — planetary 
systems: protoplanetary disks — stars: pre-main sequence 
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1. INTRODUCTION 



On the road from molecular clouds to planetary sys- 
tems, transitional disks are among the brightest sign- 
posts. Encircling T Tauri and Herbig Ae/Be stars hav- 
ing ages of 1-10 Myr, these disks have large inner holes 



tributions (SED 
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2009) and ima 
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Ratzka 
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Hug 


res et al.||2007 


Brown et al.|2009 


), tran- 



Transitional disks are so named ( Strom et al. |1990 ) be- 
cause they might represent an evolutionar y link between 
optic ally thick disks without holes (e.g., Watson et al. 
2007 1 and debr is disks conta ining only rings ot optically 
thin dust (e.g 



Wyatt 



2008). They are of special inter- 



est not least because their central clearings may harbor 
nascent planets, potentially detectable against relatively 
weak backgrounds. 

1.1. The Need for Companions 

The idea that transitional disk holes are swept clean 
by companions, possibly of planetary mass, is natural. 
We adhere to this interpretation, although not all our 
arguments as given below are the ones usually discussed. 

1.1.1. Stellar-mass Companions 

Roughly hal f of all transition al disks — 6 out of 13 in 
the sample of Kim et al. ( 2009 ) — are already known to 



contain stellar-mass companions. A prototypical exam- 
ple is CoKu Tau/4: a transition al system whose hole is 



practically empty, both of dust (D'Alessio et al. [20051 
and gas (G. Blake, private communication, 2007). Inside 
its hole of radius ~10 AU resides a nearly equal mas s K- 



the b inary can easi 



star binary haying a projected separation of 8 AU (Ire- 
land & Kraus 2008). Gravitational torques exerted by 
y counteract viscous torques in the 
staving off accretion 



disk (|Goldreich fc Tre mainc 1980jL_ 
onto either star (| Artymowicz fc Lubow 
20051. Indeed stellar accretion rates 



1994 

or CoKu 



are u nmeasurably small, < 10 M Q yr (Najita et al. 



2007) 



Ochi et al. 
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Not every hole, however, is as empty as that of CoKu 
Tau/4. In many cases there is a sprinkling of dust: opti- 
cal depths at 1 ^im wavelength for man y sources range 
from 0.01-0.1 ( | Calvet et aL||2002l |2005| . Observations 
of rovibrational emission from warm, optically thick CO 
imply that gas fills many disk holes ( Salyk et al.| [2007), 
albeit with surface densities that may be tar below those 
of conventional disks (we will argue below that this is in 
fact the case) . Most germane to our work, the host stars 
of many transitional systems actively accrete, at rates 
that on average are somewhat low er than those of con- 
ventional disks ( Najita et al.||2007[ ), but which in several 
instances approach lO -s M0yr -i , the median T Tauri 
rate. The holes of these systems must contain accreting 
gas. 

How can we reconcile the fact that many holes contain 
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gas, accreting at rates approaching those of conventional 
disks, with the fact that the holes contain only trace 
amounts of dust? To explain the paradox of simultaneous 
accretion and hole transparency, appeals arc sometimes 
made to grain growth, or the filtering of dust out of gas 



by hydrodynamic mech anisms ( Paardekooper fc Mellema 
2006[ [Rice et al. | |2006| ) or radiation pressure ( |Chiang fc: 
Murray-Clay||2007| )~ 



These proposals, which invoke changes in disk opacity, 
may be part of the solu tion. But the y cannot alone ex- 
plain the observations. Ward (2009) has criticized the 
hydrodynamic filter. The force of radiation pressure de- 
pends on uncertain optical constants and grain porosi- 
ties, and is li kely to expel only grains having a narrow 
range of sizes ( Burns et al.|1979 ) . Even if grains have the 
right properties to be blown out by radiation pressure in 
vacuum, the inward flow of accreting gas may be strong 
enough to carry as much as half of the grains that lea k 
from the rim into the hole ( |Chiang fc Murray-Clay|2007| ). 
Grain growth does not explain why transitional disk ac- 
cretion rates M tend to be severa l times smaller than 
for conventional T Tauri systems ( |Najita et al. 2007). 
Finally, none of these proposals predicts gapped ("pre- 
transitional" ) disks, which are optically thic k at stel- 



locentric dis tances a < 0.15 AU (e.g., LkCa 15; Espaillat 



et al. 2007a). 



An alternative explanation for why holes can be si- 
multaneously transparent and still contain accreting gas 
involves the special way in which disk gas accretes in the 
presence of companions, particularly those on eccentric 
orbits. If the hole rim leaks gas — and for some combina- 
tions of disk viscosity, disk press ure, and binary param- 
eters the rim can be quite leaky ( Artymowicz fc Lubow| 
[1996] |Ochi et al. 12005] — the gas can suddenly plunge in- 
ward at rates approaching freefall velocities. The catas- 
trophic loss of angular momentum is enabled by the non- 
axisymmetric and time-dependent potential of the eccen- 
tric binary, which directs gas stre amlines onto radial or- 
bits th at may intersect and shock. |Artymowicz fc Lubow| 
( |1996[ ) explained the paradox of simultaneous accretion 
and transparency: 

"... [for some circumbinary disk parame- 
ters] gravitational resonant torques are able 
to open a fairly wide gap [hole], while con- 
currently the accretion flow proceeds through 
that gap in the form of time-dependent, well- 
developed or efficient gas stream(s) carrying 
virtually all the unimpeded mass flux. The 
radial velocity of the stream is of order [the 
Kepler velocity], i.e., ~i?e [Reynolds number] 
times faster than in the disk. By mass conser- 
vation, the axially averaged surface density 
must differ by a factor of Re > 10 3 between 
the gap and the disk edge region. . . . The 
spectroscopic ramification of this is a deficit 
of the observed radiation flux emitted at tem- 
peratures appropriate for the gap location." 
(italics theirs) 

Thus reductions in the dust-to-gas ratio by grain 
growth or dust filtration are not the only processes that 
can render accreting gas transparent in transitional disk 
holes. Companions can accelerate disk gas to such high 



radial speeds that, by mass continuity, the surface den- 
sity in both gas and dust is reduced by orders of magni- 
tude. According to this explanation, the reduction in to- 
tal surface density is not necessarily due to consumption 
of gas by companions, but is rather due to gravitational 
forcing. 

1.1.2. Multiple Planetary Mass Companions 

Though a stellar-mass companion can exert torques 
strong enough to maintain holes of large size, a single 
planet-mass companion on a circular orbit cannot do the 
same job. For observationally reasonable values of the 
disk viscosity, a single companion having of order ~1 
Jupiter mass embedded within a disk at least a fe w times 



more massive carves out only a narrow gap (e .g., Lubow 



fc D'Angelo|[2006| |Crida fc Morbidellil|2007p . Disk gas 
tunn els past the planet by traveling on horseshoe-like or- 
bits ( Lubow et al.| [l999). A single planet may siphon off 
some of the gas that flows past it, but the disk accretion 
rate inside the planet's orbit is reduced from tha t out- 
side by a factor of < 10 ( |Lubow fc D'Angelo|2006"l ). This 
modest reduction in M, combined with the narrowness 
of the gap seen in simulations (Ar/r ~ 0.1), implies in- 
ner disks far too extensive and optically thick to explain 
transitional systems. 

How, then, do planets fit in? In those cases where 
the central stars of transitional disks lack stellar-mass 
companions, the same tasks of maintaining hole rims and 
increasing accretion velocities v (but not accretion rates 

M) can be performed, not by a single planet, but by 
a system of multiple planets. We imagine a series of 
planets, with the outermost lying just interior to and 
shepherding the hole rim. Gas that leaks from the rim 
is torqued from planet to planet, all the way down to 
the central star, its optical depth decreasing inversely as 
its radial speed. The more massive the planets and the 
more eccentric their orbits, the fewer of them should be 
required. 

Such a picture is supported by numerical simulations 
of Jupiter and Saturn embedded within a viscous disk 
( |Masset fc Snellgrove|200T) IMorbidelli fc Crida|2007l ). In 



these simulations the two planets were close enough that 
their gaps overlapped. Gas outside Saturn's orbit exe- 
cuted half a horseshoe turn relative to Saturn, and then 
another half-horseshoe turn relative to Jupiter, thereby 
crossing from the outer disk throug h the Jupiter-Saturn 
comm on gap into the inner disk. [Morbidcll i fc Crida| 
(20071 found that the surface density in the gap region 



was reduced by 1-2 orders of magnitude, at least near 
Jupiter. 

As our paper was being reviewed, we became aware 
of planet-disk simulations by Zhu et al. (2010, submit- 
ted) which included as many as 4 J upiter-mass planets 
and w hose results supported those of Morbidc lli fc Crida| 
( 2007 1 . Depending on the assumed efficiency with which 
planets consumed disk gas, a set of four planets was 
found to reduce surface densities in their vicinity by up 
to 2 orders of magnitude, while disk accretion rates were 
reduced by factors < 10 (see their run P4A10). However, 
such surface density suppressions are not by themselves 
large enough to explain the observed low optical depths 
of disk holes. Zhu et al. (2010) concluded that reduc- 
tions in gas opacity by some means of dust depletion 
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(e.g., grain growth) are still required. 

Companions can also accommodate gapped or "pre- 
transitional" disks in which optically thin holes contain 
optically thick annuli. As inferred from spatially unre- 
solved spectra, these annuli are narrow and abut their 
host stars, extending m ere fractions of an AU in radius 



bital frequency) 



(Espaillat et al. 2007a but see also Eisner et al. 2009 



who showed using spatially resolved observations that 
the gapped disk interpretation of SR 21 is incorrect). In 
regions far removed from secondary companions — in par- 
ticular, in those regions closest to the primary star where 
the potential is practically that of a point mass — the in- 
fall speeds of accreting gas must slow back down to the 
normal rate set by disk viscosity. By continuity, the sur- 
face density must rise back up, and optical thickness is 
thus restored. 



1.2. Companions are Not Enough: The Case for the 
Magnetorotational Instability for the Origin of Disk 
Viscosity 

Our case for companions presumes a source of disk 
viscosity. While a stellar-mass companion or a system 
of multiple planets can transport gas quickly, effectively 
generating an enormous viscosity in their vicinity (i.e., 
inside the hole), they cannot cause the outer disk to dif- 
fuse in the first place. An inviscid outer disk will not 
leak. Another source of viscosity has to act in the outer 
disk, causing it bleed inward and supply the observed 
accretion rates M . We now turn to the main subject 
of this paper, the possibility that the magnetorotational 
instability is the source of viscosity in the outer disk. 

The magnetorotational instability (MRI) amplifies 
magnetic fields in outwardly shearing disks and drives 
turbulence whose Maxwell stresses transport angula r mo- 
mentum o utward and mass inward (for a review, see Bal- 
|bus||2009|). Gas must be sufficiently well ionized for the 



MRl to operate. For the most part, T Tauri and Her 
big Ae disks are too cold at their midplanes for ther- 
mal ionization to play a role there. The hope instead is 
that X-rays emitted by host stars can provide the requi- 
site ionizati on in irradiated disk surface layers ( Glassgold 
et al.|[l997[). The basic picture was conceived by ~~ 



Gam- 



mie| (|l996), who proposed that disk surface layers ionized 
by some non-thermal means may accrete, leaving behind 
magn etically "dead" midplan e gas. Like other workers 
(e.g., Bai fc Goodman |2009 hereafter BG; and Turner 
et al. 2010[ hereafter ICS), we focus in this study on 
ionization of H 2 by X-rays. Ionization of trace species by 
ultraviolet (UV) radiation is also potentially important — 
we discuss this topic briefly at the close of our paper. 

The exposed rim of a transitional disk constitutes a 
kind of surface layer. X-rays may penetrate the rim wall, 
activate the MRI there, and disl odge a certain radial col- 
umn of gas every diffusion time ( |Chiang &: Murray-Clay 



2007| hereafter CMC). Within the MRl-active column, 
both the magnetic Reynolds number 



Re 



~D 



10 



-13 



T 



100 K 



1/2 



VlO- 8 / Vl0 10 cm- 3 / \AUJ 



3/2 



Am ■ 

u v in - / v idi" ,■!>>■■■•' / \ \n / 

( 2 ) 

must be sufficiently large for magnetic fields to cou- 
ple well to the overwhelmingly neutral disk gas. Here 
T is the gas temperature, c s is the gas sound speed, 
h = Cs/fi is the gas scale height, fl is the Kepler or- 
bital frequency, D = 234 (T/K) 1 / 2 x^ 1 cm 2 s" 1 is the 
magnetic diffusivity, x e rn is the fractional abundance of 
electrons (ions) by number, nn 2 is the number density 
of hydrogen molecules, /3- m s» 1.9 x 10~ 9 cm 3 s _1 is the 
collisional rate coef ficient for ions to s hare their momen- 
tum with neutrals ( Draine et al.||1983 ), and a is the disk 
radius. 

Dimensionless number ([I]) governs how well magnetic 
fields couple to plasma, while (pi assesses how well 
plasma couples to neutral gas. Both these numbers must 
be large for good coupling between magnetic fields and 
neutral gas. Numerical simul ations have suggested crit- 
ical values Re* of ~10 2 -10 4 ( |Fleming et al.||2000| p] de- 
ending on the initial field geometry, and Am* of ^10 2 



Hawley k, Stone|1998 hereafter HS). Some studies (e.g. 
TC5J assumed Am* ~ 1 in their determination of the 
thicknesses of MRI-active surface layers, but numerical 
simulations of marginally coupled ion-neutral systems in- 
dicated Am* may be 2 orders of magnitude higher (HS). 
The value of Am* is critical to our work. 

For typical T Tauri parameters, CMC found active ra- 
dial column densities N* ~ 5x 10 23 cm -2 or equivalently 
mass columns of S* ~ 2 g cm -2 — essentially the stopping 
column for 3 keV X-rays. When they combined their de- 
rived value for N* with an assumed value for the dimen- 
sionless disk viscosity a ~ 10 -2 , the accretion rates of 
many transitional systems were successfully reproduced. 
According to this model, the maximum accretion rate 
M inside the hole is set by conditions at the rim wall, 
i.e., by how large a radial column N* the MRI can draw 
from the rim. Stellar or planetary comp anions, known 
or suspected to be present (Section 1.1 1, regulate how 
quickly this leaked material spirals onto the host star — 
these companions modulate the radial inflow speed v{a) 
and thus the surface density £(a) = M / (2irva). But the 
companions inside the hole do not initiate disk accretion. 
They may reduce M by exerting repulsive torques to keep 
material in the rim wall from leaking in, or by accreting 
material that flow s past (e.g., Lubow & D'Angelo 2006 



Najita et al.|2007 ). But they do not generate a non-zero 



and the ion-neutral collision rate (normalized to the or- 



M in the first place. That fundamental task is left to the 
MRI operating at the rim — or whatever source of anoma- 
lous viscosity must be present in the outer disk to make 
it bleed. 

1 The magnetic Reynolds number as we define it is not 
as accurate a predictor of MRI turbulence as the Elsasser 
(a.k.a. Lundquist) number, which is given by fTJ) with c s replaced 
by the vertical Alfven speed n a y ■ Self-consistent resistive MHD 
simulations b y|Turner et al.] ( |2007^ found that MRI-active r egions 
coincide with Kls asscr numbers g rea ter than unity (see also |Sano| 
|fe Inutsuka|200l||Sano fc Stone|2002| . Our criterion Re > 10^-10* 
offsets some of the inaccuracy because v\ z < 10 c s in simulations 
of MRI turbulence. In any case we will find that the limiting factor 
for active surface layers is not Re but rather Am. 
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1.3. The Threat Posed by Polycyclic Aromatic 
Hydrocarbons to the MRI 

One concern raised by CMC but left quantita- 
tively unaddressed is the degree to which ultra-small 
condensates — macromolecules whose sizes are measured 
in angstroms — may thwart the MRI. In planetary atmo 
spheres, aero sols can strongly damp electrical conduc 



tiviti es (e.g., Schunk fc Nagy||2004 



Borucki & Whitten 



2008 )F] Most studies of active layers neglect aerosols 
and fixate on roughly micron-sized grains, despite the 
fact that in many particle size distributions, the small- 
est particles collectively present the greatest geometric 
surface area and therefore the greatest cross section for 
electron adsorption and i on recombination. Exceptions 



include Sano et al. (2000), who in one model considered 
a grain size distribution extending down to 0.005 /im = 
50 A, and BG, who considered grain sizes as small as 
0.01 /im = 100 A. Both studies found that in principle 
small grains can be deadly to the MRI. 

Notwithstanding their possibly decisive role, small 
grains are sometimes wishfully dismissed as being de- 
pleted in number by grain growth, i.e., assimilated into 
larger grains. Undeniably gr ains grow (Blum & Wurm 
|2008[ |Chiang fc Youdin|2010 1 , so much so that their col- 
lective mass may be concentrated in particles millimeters 
in size. But the question relevant for ionization chemistry 
is not where the mass is weighted in the size spectrum of 
particles, but rather where the collective surface area for 
charge neutralization is weighted. Determining the grain 
size distribution in disks seems a problem that canno t be 
forward modeled with confidence. Sano et al. ( |2000[ ) and 
BG instead parameterized the population of small grains 
and studied the effects of varying their numbers, leav- 
ing undecided the question of whether their parameter 
choices were favored by ob servation or theory. 



Like Sano et al. 



( |2600[ ) and BG, this paper consid 
ers the effects ot small condensates on the MRI. What 
is new about our contribution is that we consider the 
smallest imaginable condensates that are still accessi- 
ble to observation: polycyclic aromatic hydrocarbons 
(PAHs). These molecules, typically containing several 
dozens of carbon atoms, are excited electronically by 
ultraviolet radiation and fluoresce vibrationally at 3.3, 
6.2, 7.7, 8.6, 11.3, and 12.7 /im, the signat ure bands of 



their constituent C-C and C-H bonds (e.g., Li & Draine 



2001 Pendleton & Allamandola 2002). Spitzer satellite 



spectra and ground- based adaptive optics imaging reveal 
PAHs to be fluorescing strongly in Herbig Ae/Be and T 
Tauri disk surface layers directly exposed to stellar ul- 



travio let radiation (Geers et al. 2006, 2007 Goto et al 



20091) . Thus PAHs help to constrain the aerosol abun- 
dance where magnetically driven accretion is thought to 
occur: in disk surface layers. 

In this work we incorporate PAHs into a simple chemi- 
cal network to assess the proposal that X-ray driven MRI 
operates in disk surface layers, either on the top and bot- 
tom faces of conventional hole-less disks, or at the rims 
of transitional disks. We make as realistic an estimate as 

2 Some fire alarms work on this principle. A radioactive source 
inside the alarm drives ionization currents in the air which normally 
complete an electrical circuit. When smoke particles from a fire 
reduce the density of free ions and electrons in air, the circuit is 
broken and the alarm is triggered. 



we can of the PAH abundance based on observations, to 
gauge how deep the X-ray-irradiated, MRI-active layer 
might actually be. 

To summarize this introduction: companions — either 
stars or a system of multiple planets, but not a single 
Jupiter-mass planet — can help clear the extensive holes 
of transitional disks. The outermost companion serves to 
establish the location of the rim where viscous torques in 
the disk and gravitational torques from the companion 
seek balance. If gas leaks inward from the outer disk, it 
is driven onto the host star so quickly by gravitational 
torques from companions that its optical depth may be 
reduced by orders of magnitude. Companions, together 
with reductions in disk opacity by grain growth, thus 
maintain the transparency of the hole while still permit- 
ting stars to accrete gas. But companions do not, in 
and of themselves, cause gas in the outer disk to dif- 
fuse inward. That responsibility may be reserved for the 
MRI — whose ability to operate despite the presence of 
charge-neutralizing PAHs is the subject of this paper. 

Our paper is organized as follows. The ingredients of 
our numerical model for X-ray-driven ionization chem- 
istry in disk surface layers are laid out in Section [2j 
There we gauge what PAH abundances in disks may be. 
Results — principally, how Am and Re vary with the col- 
umn density penetrated by X-rays, and the extent to 
which PAHs reduce these numbers — are presented in Sec- 
tion [3j Analytic interpretations of our numerical results, 
and direct comparison with previous calculations (BG, 
TCS), are given there as well. We discuss our main re- 
sults for X-ray driven MRI in Section |4j and close by 
discussing the possibility of UV-driven MRI. 

2. MODEL FOR DISK IONIZATION 

In this paper we are interested in the degree to which 
stellar X-rays and Galactic cosmic-rays can ionize H2 gas 
in T Tauri disks. In this respect our study is similar 
to many others, and we make di rect comparisons of our 

For simplicity our 



work to BG and TCS in Section 3.4 



model neglects ionization of trace species like C and S 
by stellar UV radiation. Omitting UV-driven chemistry 
renders our model inconsistent because our model also 
incl udes P AHs, whose abundances we constrain in Sec- 
tion 2.4.3 by using observed PAH emission lines excited 
by stellar UV radiation. We wil l discuss the critical issue 
of UV ionization in Section 14.11 




Fig. 1. — X-ray ionized surface layers, located either on the top 
and bottom faces of a flared disk, or at its inner rim. Our calcu- 
lations apply to both situations, although they are more accurate 
for the former. We assume in this work that the lengthscale I over 
which gas is distributed radially at the hole rim is equal to h, the 
vertical gas scale height. Other sources of ionization are interstel- 
lar cosmic-rays and ultraviolet radiation from the star. At large 
stellocentric distances (a > 30 AU), cosmic rays may penetrate the 
disk from the side. For most of our paper, we neglect ion izat ion 
of trace species by far ultraviolet radiation, but in Section | we 
briefly discuss this important topic. 
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2.1. X-ray and Cosmic-ray Ionization Rates and Gas 
Densities 

Chandra spectra of pre-main-sequence stars in the 
Orion Nebula can be fitted by a pair of thermal plas- 
mas with characteristic temperatures kTx ~ 1 and 3 
keV, where k is Boltzmann's constant, a nd compara- 
ble luminosities L x ~ 10 28 -10 31 erg s" 1 (|Wolk et al. 



magnitude higher than that of the solar wind today, seem 
likely to shield disk surfaces fr om cosmic-rays directe d 



2005 Preibisch et al. 2005). The softer component is 
believed" to be emitted by shock-heated accreting gas 
( Stelzer fc Schmitt|2004[), and the harder com ponent by a 



stron gly magnetized and active stellar corona ( Wolk et al. 
2005 1 . X-ray luminosities tend to increase with increas 



mg stell ar mass and decreas ing accretion rate (see Fig- 
ure 1 of|Telleschi et al.|2"007l and Figure 17 of |Preibisch| 



et al. 2005p . These correlations are relevant for transi- 
tional disks because some transitional disks are hosted 
by higher mass Herbig Ae stars, and accretion rates for 
transi tional disks tend to be lower than for conventional 
disks (iNajita et al. 2007). For our standard model we 
lf? y 



will adopt Lx 



but we will also experi- 



ment with Lx = 10 31 erg s 1 . We fix the temperature 



erg s 
erg s" 

of the X-ray emitting plasma at kTx — 3 keV, an as- 
sumption that ignores how X-ray s pectra harden with 
increasing Lx ( Preibisch et al.||2005[ ). Although none of 
our numerical models explicitly considers kTx > 3 keV, 
we will discuss quantitatively in Section [3.3.1 how our re- 
sults scale with the higher ionization rates afforded by a 
harder (kTx — 8 keV) X-ray spectrum; we will see there 
that the effects are not largej^j See Table [l] for a list of 
all model parameters. 
We derive X-ray ionization rates Cx as a func t ion of 



penetration column N from Igea & Glassgold ( 1999 
hereafter IG), who constructed a Monte Carlo radiative 
transfer model that accounts for Compton scattering and 
photoionization. Compton scattering enables X-ray pho- 
tons to penetrate to deeper columns than would other- 
wise be possible. For our standard model we use IG's 
Figure 3 for a thermal plasma of Lx = 10 29 erg s _1 and 
kTx — 3 keV. Their ionization rates were computed for 
stellocentric distances of 5 and 10 AU; we scale these 
rates to the stellocentric distances of our model, a = 3 
and 30 AU, using the geometric dilution factor aT 2 . We 
test the accuracy of this approach by scaling their results 
internally using this dilution factor, finding (by neces- 
sity) excellent agreement at low columns where material 
is optically thin to X-rays, and agreement better than a 
factor of two at the highest columns calculated. For the 
case Lx = 10 31 erg s _1 , we increase all ionization rates 
from our standard values by a factor of 100. 

Interstellar cosmic-rays can also ionize disk gas, but are 
attenuated by magnetized stellar winds blowing across 
disk surface layers. Even the contemporary solar wind, 
characterized by a mass loss rate of ^1O _14 M yr , 
modulates the cosmic-r ay flux at Earth by as muc h as 
~10% with solar cycle ( |Marsh fc Svensmark||2003| ). T 
Tauri winds, having mass loss rates up to 5 orders of 

3 A minority of sources surveyed by Chandra exhibited superhot 
X-ray flares with peak Ly ~ 10 32 erg s - 1 and kTx 15 keV 
| |Getman et al.|2008a] |Getman et al.|2008b) . The degree to which 
superhot tlares enhance ionization rates depends on the uncertain 
flare duty cycle. Because only <~10% of the Chandra sources flared 
once or twice over a 15-day observing period, and because each 
flare lasted less than ~1 day, the extra ionization from superhot 
flares may be modest. 



normal to the disk plane (cf. Turner & Drake 2009) 



Nevertheless, cosmic-rays may reach disk gas from the 
"side," striking the disk edge-on from the outside. At 
a = 3 AU we estimate that these "sideways cosmic-rays" 
are too strongly attenuated by intervening disk gas to be 
significant. The same is not true on the outskirts of the 
disk at a = 30 AU, where column densities measured ra- 
dially outward may be sm aller than the cosmic-ray stop - 
ping column of 96 g cm" " 2 ( |Umebayashi fc Nakano|l98l [). 
Thus we consider another model at a = 30 AU where in 
addition to our standard X-ray source we include side- 
ways cosmic-rays with a constant, column -independent 



ioniza tion rate of Ccr ~ (1/4) x 10 s (Casern et al 



19981 hereafter C98) . The factor of 1 /4 is approximately 
the traction of the celestial sphere (centered on disk gas 
at 30 AU) that is not shielded by stellar winds. The total 
ionization rate C = Cx + Ccr- 

In all our simulations we neglect ionization b y energetic 
protons emit ted by the stars. As discussed by |Turner fc 
|Drake (2009), estimates of the stellar proton flux rely on 
extrapolated scaling relations, and the ability of parti- 
cles to reach the disk surface in the face of strong stel- 
lar magnetic fields is uncertain. Moreover, protons are 
emitted in flares which may occur too infre quently to 
sustain disk ionization. In one of their models, | Turner fc 
|Drake| ( |2009[ ) used a time-steady stellar particle luminos- 
ity whose ionization rate exceeded, by a factor of 40 at a 
mass column of £ = 8 g cm -2 , that of an X-ray source 
having L x = 2 x 10 30 erg s" 1 and kT x = 5 keV. This 
model probably yields a hard upper limit on the stellar 
proton ionization rate, derived under a set of generous 
assumptions. Our Lx — 10 31 erg s _1 case produces ion- 
ization rates C that approach those of the aforementioned 
model to within an order of magnitude. In any case we 
will see in Section 13.3.11 how our results can be scaled to 
anyC- 

Figure [l] depicts schematically how X-rays irradiate 
disk surface layers, usually pictured in the vertical di- 
rection as ensheathing the disk on its top and bottom 
faces. But in a transitional disk, a surface layer may also 
be present in the radial direction, at the rim of the central 
hole. We consider each of these environments in turn, es- 
timating local number densities uh 2 [Efe cm -3 ] from the 
column density N [H2 cm -2 ] penetrated by X-raysj^] 

2.1.1. Surface Layers I: Top and Bottom Faces of a 
Conventional Flared Disk 

When considering the surface layers of a conven- 
tional non-transitional disk, we describe our results as 
a function of the vertical column density N of hydro- 
gen molecules, measured perpendicular to and toward 
the disk midplane. Thus our N coincides with N± of 
IG, save for a factor of 2 because IG count hydrogen nu- 
clei whereas we count hydrogen molecules. An equivalent 
measure of vertical column density N is the mass surface 
density £ = Nfi, where /1 w 4 x 10~ 24 g is the mean 
molecular weight of gas. 

4 In this paper, ionization rates f , column densities N, and frac- 
tional densities x are referred to hydrogen molecules, not hydrogen 
nuclei. 
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TABLE 1 
Model parameters. 



Parameter 



Variable 



Value 



Reference 



Disk radius a 

X-ray source luminosity" Lx 

X-ray source temperature &Tx 

Cosmic-ray ionization rate Ccr 

Initial CO abundance 6 %co 

Total metal abundance* 1 ' 6 xm 

Total grain abundance 6 ^grain 
Grain settling (depletion) factor e gr ain 

Total PAH abundance 6 £pAH 

PAH depletion factor e PAH 

Central stellar mass M* 

Gas temperature T 



3, 30 AU 

10 29 (10 31 ) erg s" 1 
3 keV 

0, (0.25 x 10" 17 ) s" 1 
10" 4 

10" 8 (0, 10~ 6 ) 
6 x 10" 15 
lO" 3 < 



grain 

in < io- 1 



c gram 

10- B epAH 
10" 5 < epAH < 
1 Mq 
80, 30 K 



2.1 



Section 
Sectio n Iz.i 
|CaselIi eTlu"] ||199 8|| 

Aikawa et al. ( 1£ 
Section 
Section 
Section 
Section 
Section 



2XT 



TF2 



Section 



a Values in parentheses correspond to test cases different from our standard model. 
b All abundances are relative to H2 by number. 



To good approximation, the local number density 

n H2 « N/h (3) 

where the vertical scale height h = c s /0 = (fcT//i) 1 / 2 /J7. 
For gas temperature T rj 80(30) K at a = 3(30) AU (sec 
Section 2.2 for how we derive these temperatures), we 
find h =TTD9(1.8) AU. 

Equation pi underpins all our calculations of chemical 
equilibrium. For a typical N ~ 10 23 H2 cm~ 2 , we have 
n ~ 7 x 10 10 (4 x 10 9 ) H 2 cm" 3 at a = 3 (30) AU. 

2.1.2. Surface Layers IT. Gap Rim of Transitional Disk 

We assume that the gap rim is not shadowed from 
the star by gas interior to the rim. We cannot prove 
that the rim is not shadowed, but disk models based on 
the i nfrared SED suggest it is not (e.g 
2005 1 . Possibly gas at the rim wall 



Calvet et al. 
puffs up" because 

it is heated by X-rays and can maintai n a larger vertical 
heigh t than gas inside the hole (e.g., Dullemond et al. 
200l| . 

For the case of the rim of a transitional disk, we reinter- 
pret N (equivalently E) as the radial column of hydrogen 
molecules traversed by X-rays (Figure 1). To estimate 
the local number density nn 2 , we need to know the ra- 
dial lengthscale £ over which materi al at the rim wall is 
distributed. Plausibly h < I < a. 
Clayl (12007} take £ ~ a, but models 



Chiang fc Murray- 
jased on SEDs and 
images suggest the rim is much sharper. Here we assume 
that £ ~ h so that Equation ^ applies equally well to 
transitional disks as to conventional disks — keeping in 
mind that N should be measured radially for the former 
and vertically for the latter. 

To calculate ionization rates in the rim, we still use 
the results of IG, reinterpreting their N± in their Figure 
3 as our radial column N. Clearly the scattering geom- 
etry differs between the case of a transitional disk rim 
and the case of the top and bottom faces of IG's con- 
ventional disk. Where material is optically thin to stel- 
lar X-rays, the two cases match in ionization rate (per 
molecule), but where it is optically thick, we underes- 
timate the ionization rate in transitional disk rims by 
using IG because more X-rays escape by scattering ver- 
tically out of conventional disk surface layers than from 
the rim. Another reason we underestimate the ionization 
rate at high column density is because the total X-ray 



flux per unit surface area of the disk is lower for conven- 
tional surface layers — which are illuminated at grazing 
incidence — than for the rim, which is illuminated at nor- 
mal incidence. Nevertheless we estimate that these er- 
rors are of the order of unity, insofar as the columns that 
might possibly be MRJ -active are not too optically thick 
to X-rays (Section 3.3 1, because the Thomson scattering 
phase function is fairly isotropic, and the cross-section 
for scattering is only comparable to that for photoion- 
ization at the relevant photon energies. In any case we 
will explore the effects of higher ionization rates by run- 
ning a model with higher Lx 
3.3.11. 



10 3 erg s 1 (Section 



2.2. Gas Temperature 

The surface layers of protoplanetary disk atmospheres 
vary widely in temperature, from ^5000 K at the lowest 
columns where stellar X-rays heat the gas, down to < 100 
K at the highest columns where dust reprocesses optical 
starlight. We draw our temper atures from the thermal 
model of Glassgold et al. ( 2004 ), which in turn i s based 
on the dust temperature model of D'Alessio et al. ( 1999 1. 



At a = 1 AU at column densities of interest (N > 10 2 ' 
cm -2 ), thermal balance is controlled primarily by repro- 
cessing of starlight by dust, and gas and dust tempera- 
tures are nearly equal at ~130 K ( jGlassgold et al.|[2004 
their Figure 2). We adjust this result for the disk radii o 
our standard model using the dust temperature scaling 
law f or the midplane of a passiv e flared disk, T cx a -3 / 7 

l|l997|). Thus at a = 3 AU we 
30 K. 



(e. 

have 



Chiang & Goldreich 



K, and at a 



"3Tj AU we have T 
Note that these temperatures are lower — and arguably 
more realistic — than those assumed by BG and TCS, who 
invoked temperatures of the traditional Hayashi nebula 
without justification. 

2.3. Chemical Network 

Following [rlgner fc Nelson] ( |2006[ hereafter IN), CMC, 
and BG, we apply a simple network of chemical r eactions 
based on that designed fo r molecular clou ds by lOppen- 



heimer fc Dalgarno (1974 hereafter OD). Ilgner & Nel- 



son| ( |200fc>| ) and BG compared the results of UU- based 



schemes to those of more complex networks extracted 
from the UMIST (Unive rsity of Manchester In stitute of 
Science and Technology; |Woodall et al.| [20071 hereafter 
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7; i Vasyunin et al. 2008) database. Fractional elec- 

. am; 



W07: 

tron abundances derived by IN using the simple network 
were greater than those derived using the complex net- 
work, whereas BG, who used a more recent version of the 
UMIST database, found that the sign of the difference 
varied from case to case. The magnitude of the difference 
ranged up to a factor of 10, but was often < 3. Using the 
simple network seems the most practical approach, if we 
are content with order-of-magnitude answers. In Section 
|3.4[ we test the results of our code against those of BG 
^Sd TCS. 

All reactions in our OD-based network are listed in 
Table [2] and shown schematically in Figure [2j Rate co- 
efficients and their temperature dependences are taken 
from the UMIST database. The chain of events basi- 
cally proceeds as follows. X-rays ionize H2 to H^, which 
rapidly reacts with H2 to produce . The Hjj" ion com- 
bines with CO to form HCO + . Most HCO + ions disso- 
ciatively recombine with free electrons, but some transfer 
their charge to gas-phase metal atoms such as Mg. Free 
metals tend to be abundant positive charge carriers, as 
they recombine with free electrons only by a slow radia- 
tive channel. Charged particles in the network (e~, Hj", 
HCO + , metal + ) can neutralize by collisionally transfer- 
ring their charge to PAHs and grains. The collisional 



charging process is described in Section 2.5 



The reaction loop is closed by the formation of H2 on 
grain surfaces. To compute the rate of this reaction, we 
take neutral H atoms to collide with grains using the ge- 
ometrical cross section for grains, and adopt from BG 
the uniform probability n = 10~ 3 for a pair of adsorbed 
hydrogen atoms to form a hydrogen molecule (see their 
Equation 27). The precise rate of this reaction is not im- 
portant for us, as it only sets the equilibrium abundance 
of H, which is irrelevant for the ionization fraction, as 
long as n H <C n H2 ■ 

In their original study OD included ionization of He 
and reactions involving atomic and molecular oxygen. 
We neglect these for simplicity. Most reactions involving 
oxygen initiate with the formation of the hydroxyl ion 
(H3" + O — > OH + ), which proceeds at a rate only com- 
parable to the formation of HCO + , which we do account 
for. Thus our neglect of oxygen within the OD frame- 
work is not expected to alter our results for the fractional 
ionization by more than a factor of 2. In any case, in Sec- 
tion |3.4| we will compare our results with those of more 
complex networks considered by BG and TCS. 

2.4. Properties and Abundances of Trace Species 

The trace in gredien ts of our model include gas-phase 
metals (Section |2.4.1 ), a monodispersion of micr on-sized 
grains (Section [5X2]), and PAHs (Section |2~4~3| ) . 



2.4.1. Gas-phase Metals (Magnesium) 
For gas-phase m etals which serve imp ortantly as elec 



tron donors (OD; Fromang et al. 2002), we are guided 
by Mg, whose sola r abundance is 3.5 x 10 -5 atoms per 
hydrogen nucleus (|Lodders||2003| ). The fraction of Mg 
that is in the gas phase — neither ncorporated into grain 
interiors nor adsorbed onto grain surfaces — might be at 
most 3-30% by number, its value in the diffuse interstel- 



lar medium (Jenkins 
of protoplanetary dis. 



2009). In the dense environments 
cs, the gas-phase fraction should be 



much smaller because magnesium is used toward building 
grains. 

Nominally, our model temperatures of 30-80 K are so 
low that almost all of the Mg not incorporated into grain 
interiors should be adsorbed onto grain sur faces, leaving 
behin d only a tiny fraction in the gas phase ( Turner et al. 
2007 their Section 2.2; see also Equation 26 of BG). Just 
how tiny is uncertain, given how sensitive the adsorption 
fraction is to gas temperature, and how stee p temper- 
ature grad ients can be in disk surface layers (Glassgold 
et al.|2004[ ). Turbulent mixing of hot, high altitude, nor- 
mally metal-rich layers with cold, low altitude, n ormally 
metal-poor layers can also complicate matters (Turner 
eFaLl[2007| TCS). 



We adopt a standard metal abundance of xm — 10 
metal atoms per H2, which corresponds to a gas-phase 
fraction of ~10 by number relative to solar. Our choice 
is similar to those of IN and BG. We also experiment 
with a metal-free case in which all metals have been ad- 
sorbed onto grain surfaces (xm = 0), and a metal-rich 
case for which im = 10 -6 . Although the metal-rich case 
is not especially realistic and is not justified by our model 
parameters — in particular our low gas temperatures — we 
consider it anyway because we would like to understand 
the effects of metals in principle, and to connect with 
other studies that consider simila rly large metal abun- 
dances (CMCjlTuTner et al.||2007| TCS). 



2.4.2. Grains 
The number of grains per H2 molecule is 

P Pdust 



^TTS 3 p s p ga 



(4) 



where fi « 4 X 10 -24 g is the mean molecular weight 
of gas and Pdust/Pgas is the dust-to-gas mass ratio. For 
simplicity we consider grains of a single radius s — 1 /im 
and internal density p s — 2gcm~ 3 . There is ample evi- 
dence that micron-sized grains abound in surface layers, 
both from mid-infrared spectra of silicate emission lines 
(e.g., Natta et al.||2007 ) and from scattered light images 
at similar wavelengths (e.g., McCabe et al.|[2003 ). 

The dust-to-gas ratio in surface layers may differ con- 
siderably from its value in the well-mixed diffuse inter- 
stellar medium (ISM): 



Pdust _ Pdust 
— Egrain 

Pgas 



(5) 



ISM 



where for the ISM o f solar abundance /Odust//Ogas|iSM = 
|2003j. Based on model fits to observed 



0.015 (Lodders 



2006 Dullemond & Dominik 2004), There is consensus 



004) 
fly 1 



far- infrared SEPs flChiang et al.||20 0l] D'Alessio et al 

nere 



that surlace layer grains directly illuminated by opti- 
cal light from their host stars have settled toward the 
midplane into regions of denser gas. Thus e gra i n < 1, 
but actual values are not known with certainty, because 
small changes in the SED resulting from small changes in 
the height of the dust photosphere imply large changes 
in gas density in a near-Gaussian atmosph ere. For ex- 



ample, Figure 15 of D'Alessio et al. (2006) shows that 
changing the far-infrared SED by less than a factor of 2 
changes eg ra i n by a factor of 10. 

Table [3] lists fitted values of e gra i n for some transitional 
disks, drawn from the literature. At best they are ac- 
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TABLE 2 

Chemical Reactions Including Collisional Charging of PAHs and Grains. 



lVrnm rifvr 

1 N Hill l-f^L 












Value 
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-> H+ + e- 
-> + c- 
-4. +H 
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radiative transfer model 


IG 


2 


H 2 + Cosmic-ray 


CCR 


(1/4) X 1Q- 


-17 g -l 


C98 


3 


H+ + H 2 


°H+ Ho 


2.1 x 10" 


-9 




W07 


4 


H+ + CO 


v HPO+ -t- Ho 


"h+,CO 


1.7x10" 


-9 




W07 
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a M,X+ 


1.0x10" 
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W07 
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W07 
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-> H 2 + H 
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" 8 (T/300K)" - 52 


W07 
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HCO+ + e- 


~> H + CO 


%CO+x 


2.4x10" 


" 7 (T/300K)" - 69 


W07 
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M+ + e- 


-» M + fe 


"M+.e 


2.8x10" 


-12 
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W07 
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Section 
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Section 
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DS 
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— ► grain(Z— 1) 


^grain.c 


Section 
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DS 
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grain(Z) + X+ 


— > grain(Z+l) 


^grain,X+ 


Section 


2.5 




DS 
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PAH(Z=-1) + PAH(Z=1) 


-» 2x PAH(Z=0) 


a PAH,PAH 


Section 


2.0 




DS 


15 


H + H + grain 


— > H 2 + grain 


^physisorption 


Section 


2.:, 




BG 



and a has units of < 



M represents a gas-phase atomic metal, e.g., Mg. 
d X+ can be cither H+ HCO+, or M+ . 

curate to order of magnitude. For our calculations we 
consider 10" 3 < e graill < 10" 1 (Table m). 



TABLE 3 

Dust Settling Parameter tp-rain FOR Some Transitional Disks 



"ST 



"HcTs 



LkCa 15 HP 3 Espaillat ct al. (2007a); Chiang ct al. (2001)" 

UX Tau A 10~ 2 Espaillat ct al. (2007a 

CS Cha 10" 2 |Espaillat efaT] |2007b 1 

GM Aur 10~ 1 |(Jalvet et aTf([2003{ 

DM Tau 10" 1 Calvet et al. (2DDSi 



2.4.3. PAHs 

For simplicity we model PAHs as spheres, each having 
a radius s — 6 A and internal density p s — 2 g cm -3 . 
Although in reality carbon ato ms in PAHs are arranged 
in sheets and not spheres (e.g., Allamandola et aLfl 999), 
the difference in cross section arising from geometry is 
only on the order of unity. Each of our model PAHs has 
about as much mass as a real PAH containing Nq = 100 
carbon atoms. A PAH of this size is estimated to be just 
large eno ugh to survive phot o-destruction around Herbig 
Ae stars fVisser et aL|2007[ ). 

The central wavelengths of PAH emission lines from 
Herbig Ae/Be (HAe/Be) disks are observed to trend 



with the eff e ctive temperatures of their host stars ( Sloan 



et al. 20051 Keller et al. 20081. This correlation mcH 



cates that PArls m disks are not merely PAHs from the 
diffuse ISM transported unadulterated into circumstellar 
environments. Rather, PAHs in disks have been photo- 
processed, their chemical bonds altered by radiation from 
host stars. Possibly PAHs are continuously created and 
destroyed by local processes, e.g., sublimation of grain 



mantles and photodestruction ( Keller et al. 2008 ) . In this 



paper we do not account explicitly tor such processes, i.e 
we do not attempt to calculate the abundance of PAHs 
from first principles. Rather we fix the abundance of 



PAHs using observations, as detailed in the remainder of 
this subsection. 

Emission from PAHs is detected in an order-unity frac- 
tion of H Ae/Be stars, but i s rarely seen in T Tauri 



stars (e.g., Geers et al. 2006 hereafter G06). In prin 



ciple this could mean that PAHs are less abundant in T 
Tauri disks, but the more likely explanation is that this 
is an observational selection effect: Herbig Ae/Be stars 
are more luminous in the ultraviolet (UV) and therefore 
cause their associated PAHs to fluoresce more strongly 
(see, e.g., Figure 9 of G06, which shows how the PAH in- 
tensity drops below the Spitzer detection threshold with 
decreasing stellar effective temperature). Another clue 
that PAHs are just as abundant in T Tauri disks as in 
HAe/Be disks is that those few T Tauri stars with pos- 
itive PAH detections tend to have unusually low mid- 
infrared continua, allowing PAH emission lines to stand 
out more clearly (G06). In other words, those T Tauri 
systems where PAHs have been detected are transitional 
systems, and their PAH abundances seem no different 
t han in their HAe/B e counterparts. 

|Geers et al. (2006) used radiative transfer models to 
fit the intensities of the 11.2 /im PAH fluorescence line 
in eight Herbig Ae and T Tauri disks, concluding that 
the PAH abundance is 10 -7 -10 -8 per H2 (see their Fig- 
ure 9). This result is highly model dependent. Perhaps 
the chief source of uncertainty lies in the grain opacity. 
Inferred PAH abundances relative to gas are sensitive 
to assumptions about the local grain size distribution 
and dust-to-gas ratio because the soft ultraviolet radia- 
tion (~1000-3000A wavelength) which causes PAHs to 
fluoresce is also absorbed by ambient grains. Thus the 
intensity of PAH emission depends on how many grains 
are competing with PAHs for the same illuminating pho- 
tons. The grain opacity, in turn, d ecreas es by orders of 
magnitude as dust settles (Section 2.4.2). Because G06 
did not account for dust sedimentation and instead as- 
sumed the dust-to-gas ratio in disks was similar to that 
of the well-mixed ISM, the PAH abundances relative to 
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H+ + e" H 2 + H 
a = 5 X 10~ 8 cm 3 s _1 



A 




grain + H + H — ¥ grain + H2 




ht + H2 —> ht + h 



a = 2 X 10" 9 cm 3 s _1 ^~ 





H 2 + hi/ -»■ + e 



+ CO HCO+ + H 2 
a = 2 X 10 -9 cm 3 s _1 



M + X + — > M+ + X 





M+ + e" M 
a = 9 X 10 -12 




^PAHs 



dMh grams 



Fig. 2. — Our chemical reaction network, derived from [Oppcnheimcr &; Dalgarno (1974). Rate coefficients as shown in this Figure are 
evaluated at T = 80 K. See Table [2] for a comprehensive list ot all modeled reactions and precise rate coefficients. 
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gas that they computed are actually upper limits. Sur- 
face layer grains have likely settled toward the midplane 
into regions of high gas density, and were the local PAH 
abundance to remain as inferred at lO^-lO -8 per H2, 
the ratio of PAH line int ensity to dust continuum would 
be larger than observed (Dullemond et ah|2007|. 



In our model the number of PAHs per M2 is 

2-PAH = £PAH X 10~ 6 , 



(6) 

where £pah < 1 measures how depleted PAHs are in 
disk surface layers relative to PAHs in the diffuse ISM 
(Li & Drainc 2001). Based on the considerations above, 
we should combine the G06 depletion factor of 0.01-0.1 
with the grain depletion factor e gra i n ~ 0.001-0.1, which 
parameterizes the reduction of dust opacity due to dust 
sedimentation. We thus estimate that 10~ 5 < £pah < 
10- 2 (Table []}. In our calculations we select the pa- 
rameter combinations (e gra ; n , £pah) = (10 _1 ,10 -2 ) and 
(egrain, cpah) = (10 -3 , 10 -5 ) which bracket the range of 
possibilities. 

A final point to consider is whether the observed PAHs 
are present at the same column depths that are rele- 
vant for X-ray driven MRI. The X-ray stopping column 
should be compared with the column that presents op- 
tical depth unity to the soft UV radiation driving PAH 
emission. In the model of G06 in which dust has not set- 
tled, photons at wavelengths of 1000-3000 A are stopped 
by submicron-sized silicate/carbonaceous grains within 
a hydrogen column of ^0.005 g cm" 2 (V. Geers, pr ivate 
communication, 2010; see also Habart et al. 2004 who 
used similar dust opacities). After we account for jrain 
settling (e gr am), the UV absorption column increases to 
~0.05-5 g cm -2 . Although model-dependent, our esti- 
mate of the UV absorption column corresponds well to 
X-ray stopping columns, and thus to columns that might 
possibly be MRI-active. 

2.5. Collisional Charging of PAHs and Grains 

Grains and PAHs are modeled as conducting spheres 
for simplicity. Electrons and ions collide with and stick 
to grains and PAHs, charging them. When the total elec- 
tron capture rate by grains and PAHs matches the total 
ion capture rate, the distribution of charges carried by 
PAHs and grains reaches dynamical equilibrium. The 
average charge state on a PAH/grain (Z) < because 
in thermal equilibrium electrons move more quickly than 
ions. Convenient and readily deri ved approximations for 



(Z) in various limits were given by Draine & Sutin (19871 
hereafter DS). We will find that (Z) ranges between — 1 
and for our PAHs of radius 6 A, while for our micron- 
sized grains (Z) m —22. The remainder of this subsec- 
tion details how we compute the electron and ion capture 
rates. 

The rates at which ions or electrons collide with PAHs 
or grains are enhanced by Coulomb focusing between 
static charges, as well as by the induced dipol e force 
( |Natanson|[l960| |Robertson fc Sternovsky] |2008[) . The 
cross sections can be derived from kinetic theory by con- 
sidering the potential between a conducting sphere of 
radius s and charge Ze, located at a distance r from a 
charge q: 



qZe 
r 



q 2 s 3 



(e.g., Jackson 1975 ). The first term is the usual monopole 
interaction, while the second arises from the induced 
dipole (image charges). For a neutral sphere, the 
velocity-dependent cross section derives from applying 
conservation of energy and momentum to the second 
term of Multiplying this cross section by either the 
electron or ion velocity, and averaging over a Maxwellian 
speed distribution at temperature T, yields the rate co- 
efficient (units of cm 3 s _1 ) 



a = ns 2 Sc 1 



2skT 



for Ze/q = (8) 



where fc is the Boltzmann constant, c = yj 8fcT / irm is 
the mean speed for either electrons of mass m = m or 
ions of mass m = mx+ , and S is the probability that the 
electron/ion sticks to the PAH/grain. We will discuss 
the sticking coefficient S shortly. 

For a charged sphere, the cross section is enhanced by 
both terms in Q. There is no analytical solution for 
this case, but DS provided the following approximate 
formulae: 

a = 7rs 2 S'c[l-Ze/((?r)](l+ v /2/(T - 2Ze/q)) for Ze/q < 

(9) 

a = tis 2 Sc[l+{AT+ZZe/ q)- 1/2 ] 2 cM-P /t) for Ze/q > 

(10) 

where r = skT/q 2 , 



1 



Ze 



qg 2g 2 (g 2 - 1) 



(11) 



(12) 



and g is the solution to the transcendental equation 

2g 2 ~l Ze 
g(g 2 -I) 2 q 

Upon colliding with a PAH or grain, the electron or ion 
sticks with probability S. For ions, we set S = S*x+ = 1 
(DS; IN; BG). For electrons colliding with PAHs, S = S c 
depends on the deta i led m olecular structure of the PAH. 



Allamandola et al. ( 1989 ) calculated how the electron 
sticking coefficient increases with both the number of 
carbon atoms and the electron affinity. The dependence 
on electron affinity is especially strong. A PAH hav- 
ing Nq — 32 and an electron affinity of 0.7 eV has 
S c «3x 10~ 5 (see their Figure 25), while the same-sized 
PAH with an electron affinity of 1 eV has S e s» 10~ 2 (sec 
page 769 of their paper). Estimated electron affinities 
of real Nc = 32 PA Hs (e.g., ovalene and hexa benzo 



2r 2 (r 2 



(7) 



coronene) exceed 1 eV. Allamandola et al. ( 1989 ) stated 
that "only for pericondensed PAHs [which are more sta- 
ble than catacondensed PAHs] containing considerably 
more than 20 C atoms will the electron sticking coeffi- 
cient approach unity." Based on these considerations, 
we take S e = 0.1 for our N c = 100 PAHs. For the much 
larger grains we set S e = 1. 

In our code, the range of charges a grain can possess 
extends from Z — —200 to +200. We have verified that 
this range is large enough to accommodate the entire 
equilibrium charge distribution, which for our /xm-sized 
grains peaks at —22 (see Figure 4| Accounting only for a 
few charges — up to \Z\ — 3 as did Sano et al. (20001, IN, 
BG, and TCS — is not necessarily adequate for micron- 
sized grains which have fairly large capacitances. For 
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PAHs we consider charges Z between —16 and +16. Most 
PAHs will turn out to have either Z = or — 1. Because 
of their smaller size, a single PAH will be less charged 
than a single grain; electrons collide less frequently with 
a negatively charged sphere as the radius of the sphere 
decreases and the Coulomb potential steepens. 

We neglect adsorption of neutral gas-phase species onto 
grain surfaces, and any mass increase of grains and PAHs 
from collisions with ions. Grain-grain and PAH-grain col- 
lisions are negligible and ignored. We do account for the 
possibility that a PAH with a single negative charge can 
neutralize by colliding with a PAH with a single positive 
charge (reaction 14 in Table [2| , though in practice this 
reaction is not significant. 

2.6. Numerical Method of Solution 

The time-dependent rate equations for the abundances 
of species are readily constructed from the reactions 
listed in Table [2j For example, the number density of 
electrons n e obeys 



dn c 



^HiC-nc^ax+.c "x+ 



x+ 



10 



-n e apAH.o "PAH 

Z=-16 
200 



Z=-200 



(13) 



where the index X + runs over reactions 7, 8, and 9 in 
Table [2j The first sum over Z occurs over the charge 
states of PAHs, while the second sum occurs over the 
charge states of grains, with rate coefficients a given in 
Section 1231 

The charge distributions of PA Hs and grain s are 



governed by recurre nce equations ( Parthasarathy 1976 
Whitten et ai~1|2007| ), e.g., for PAHs: 



drip ah, z 
dt 



: ("PAH "PAH,c) z+1 Uc 
+ £ ( npAH a PAH,X+)z_l "X+ 

x+ 

- ("PAH OtpAK,e) z n e 



£ 

x+ 



"PAH «PAH : X+ ) y "X+ 



(14) 



The right-hand side of Equation ( 14 ) accounts for all 
the ways in which PAHs of charge 2Tcan be created or 
destroyed by collisions with electrons and ions (react ions 
10 and 11 in Table [2). When Z = ±1, Equation (p) 
is supplemented by an extra loss term accounting for 
reaction 14. 

All rate equations are discretized to first order and ad- 
vanced simultaneously using a forward Eulcr algorithm 
with a fixed timestep At < lxlO" 3 s. At t = 0, all PAHs 
and grains have Z = and all hydrogen is in the form 
of H2 . In principle we could simply advance the network 
forward until the system equilibrates, i.e., until the time 
rates of change of the abundances fall below some speci- 
fied tolerance. However the reaction rates in our network 



span almost 5 orders of magnitude. Thus, our equations 
are stiff and a brute-force integration would require an 
inordinate number of timesteps. Metals are typically the 
slowest constituent to reach equilibrium because they re- 
act with electrons only slowly by radiative recombination 
(reaction 9). 

To circumvent the bottleneck posed by metals, we pro- 
ceed as follows. We run R versions of the code having R 
evenly spaced initial abundances for charged metals "m+ • 
For each run, we initially set n c = "m+ to ensure charge 
neutrality. We run each code until the abundances of all 
species drift only because of slow changes in n M + . We 
evaluate dn M + / dt at the end of each run. The equilib- 
rium value of n M + is bracketed by the two runs having 
opposing signs for dny^+ j dt. We then start a new itera- 
tion with R runs having initial metal abundances evenly 
spaced between the two bounding runs of the previous it- 
eration. In this way we refine our initial guesses for "m+ 
until we arrive at two sets of initial conditions that differ 
by less than 30%. The equilibrium value of Um+ we re- 
port lies at the intersection of the two curves for n M +(i), 
linearly extrapolated forward in time. Other variables 
(n e ,riHco+7 an d "h + ) are a ^ so extrapolated. The num- 
ber of runs R at each iteration varies from 2 to 5. 

We use the time t cq at which the two extrapolated 
curves for n M + intersect as an estimator of the equili- 
bration time of the chemical network. For t cq to be a 
robust estimator, it should be independent of initial con- 
ditions. We found that the value of t eq remained constant 
to within a factor of 3 when initial conditions varied over 
2 orders of magnitude 
pared to dynamical timescales D, 



Our values for t oq will be com- 



in Section 3.3.2 



Errors are estimated by monitoring conservation of 
charge and conservation of the total number density of 
PAHs + grains. Over 10 10 timesteps, the charge re- 
mains constant (at zero) to better than one part in 10 s , 
with similar results for the number density of PAHs + 
grains. As a test of our code, we reproduced th e normal- 
ized charge dist ribution on PAHs computed by | Jensen fc 
Thomas (1991 see their Figure 2a). 



3. RESULTS 



In Section 3.1 we describe how charges distribute 
themselves on P AHs and grains in dynamical equilib- 
rium. In Section |3.2[ we explore how the free electron 
and ion abundances vary with increasing PAH abun- 
dance. In Section |3.3[ we show what all this implies for 
the degree of magnetic coupling in disk surface layers, in- 
terpreting our numerical results whenever possible with 
simple analytic estimates. In that section we also com- 
pute timescales for the chemical network to equilibrate, 
and compare to the dynamical timescales over which the 
MRI may act. In Section [3~4} we test the validity of our 
simple network/code by seeing how closely we can repro- 
duce the results of more complex networks/codes by BG 
and TCS. 

3.1. Charge Distributions on PAHs and Grains 

Figures [3] and [4] show the charge distributions on 
PAHs and grains, respectively, for the case a = 3 AU, 
£ = 0.3gcm~ 2 , xm — 10~ 8 (standard metal abundance), 



epAH = 10 (low PAH abundance), and e g 



10" 



(low grain abundance). Most of the PAHs either have 
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Zpah = or Zpah = — 1. For grains, the average 
charge state (the peak of the distribution) is (Z gla i n ) w 
—22. The shape of the charge distribution for grains ap- 
proaches the Gaussian given by Equation (4.15) of DS. 

We may understand (Z) simply. Consider the PAHs; 
identical considerations apply to grains. We take the 
limit that the dominant ions are metals and the limit 
that the total charge carried by PAHs is much less than 
the free charge. Together these limits imply that x c f=s 
xm+ ■ Then detailed balance between forward and reverse 



reaction rates dictates that (cf . Equation 14 ) 

("PAH QPAH,c) z+1 = (™PAH"PAH,M+) Z 



(15) 



From this equation it is evident that if ever the rate co- 
efficients (apAH,e)z+i and («pah,m+)z were to be equal, 
the densities (tipah)z+i and (npAij)z would be equal, 
i.e., the charge distribution would be at an extremum. 
Thus we may estimate the average charge (Z) by merely 
plotting the rate coefficients o:pah,m+ and ^PAH.e against 
Z and seeing where the curves intersect. This exercise 
is performed in Figures [3] and [4j Indeed what the full 
numerical model gives for (Z) is close to the Z for which 
the curves for the rate coefficients intersect. (Of course, 
perfect agreement cannot be obtained because it is never 
strictly true that (a PA H,o)z+i = (o:pah,m+)z-) 

There is another, even simpler limit where (Z) may be 
estimated. In the extreme case that the gas is so satu- 
rated with grains or PAHs that practicallyno free charges 
are left, we must have (Z) — » 0. Figure [5] shows the re- 
sults of an experiment using our full code m which we in- 
crease epAH until this regime is reached. For this Figure, 
the grain abundance is set to zero to isolate the effects of 
PAHs. Figure [6] is analogous; e gra i n is increased while the 
PAH abundance is held fixed at zero. Both figures follow 
the transition from (Z) ^ to (Z) — > 0. Observationally 
inferred values for e gra ; n (see the shaded region of Figure 
6| are never so high as to cross into the (^ gra in) — > 
regime. By contrast, Figure [5] shows that PAHs may be 
sufficiently abundant in disks that they impact the den- 
sity of free charges. The critical PAH abundance x PAH 
dividing the (Zpah) limit from the (Zpah) — > limit 
is the one for which an electron attaches itself to a PAH 
as frequently as it recombines with an io n. T his critical 
abundance is discussed further in Section [3.21 



3.2. Ionization Fraction vs. PAH Abundance 

Figure [7] plots the fractional electron and ion densi- 
ties, x e and Xi, against the PAH abundance xpah, for 
a = 3 AU, £ = 0.3 g cm~ 2 , im = 10~ 8 (standard metal 



abundance), and e g 



io- 



-10 (see the figure car 



tion for how e gra i n is assigned to each €pah)- Figure I 
is identical except that it considers the metal-rich case 
xm = 10 -6 . The primary ions in both cases are atomic 
metals and HCO + molecules. At low PAH abundances, 
charged metal ions are the most abundant. As the num- 
ber of PAHs is increased, HCO + becomes the dominant 
ion. See Table [4] for a precise breakdown of component 
ion densities for our standard metal abundance case. 

According to Figures [7] and [8j the electron and ion den- 
sities are nearly equal and constant with xpah as long 
as £pah is not too large. In going from the standard 
metal abundance of Xm = 10 -8 to the metal-rich case 
of xm — 10 -6 , the free charge abundance increases by 
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Fig. 3. — Equilibrium charge distribution on PAHs (solid cir- 
cles, left axis) for a = 3 AU, S = 0.3 gem -2 , im = 10 -8 (stan- 
dard metal abundance), epAH = 10 — 5 (low PAH abundance), and 
Egrain = 10~ 3 (low grain abundance). The distribution peaks at 
Z = 0, approximately where the attachment coefficients (dashed 
lines, right axis) for electrons with PAHs and metal ions with PAHs 
cross. 
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Fig. 4. — Same as Figure [3] but for grains. 

an order of magnitude. Once xpah exceeds some critical 
abundance x PAil , the electron and ion densities diverge — 
the ion density is higher, and the balance of negative 
charges is carried by PAHs. In the limit xpah 3> ^pahj 
both the electron and ion densities decrease with increas- 
ing PAH abundance in an approximately inverse linear 
way. All of this behavior can be understood analytically 
as follows. 

3.2.1. Analytical Model for Ionization Fraction vs. PAH 
Abundance 

Our code's results for x e (xpAti), ^i(^PAH), 

and ^pah 

may be understood using the following simple model. 
The model consists only of X-rays, molecular hydrogen, 
electrons, PAHs, and one ion species — either HCO + for 
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TABLE 4 

Densities of Charged Species for our Standard Model (x m = 10 _8 ,Lx = 10 29 ergs _1 ) 
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our standard metal abundance case, or ionized metals 
M + for the metal-rich case. In the simplified model, X- 
ray ionization of a hydrogen molecule produces a free 
electron and — skipping the entire reaction chain — one 
ion. The system reduces to the rate equations 



dn c 
~dt 
drii 
~dt 



2 - n c riia iyC - n c a;pAH"-H 2 (,apAH,e)<z) 



(16) 



C™h 2 - nin c a^ c - nio;pAH«H 2 (apAH,i)(Z) (17) 



and 



The sub- 



for the electron and ion densities, n c 
script i denotes either HCO + or M + . 

In the simplified model, all PAHs with abundance 
^pah"-h 2 are assumed to be identically charged. We set 
this common charge equal to the average char ge state 
(Z), results for which were given in Section [3~T| 

We exclude our large, micron-sized grains from the an- 
alytic model. Although these grains were useful for infer- 



ring PAH abundances from observations (Section 2.4.3), 
their collective surface area is too low to significantly in- 
fluence the electron chemistry in any of our model runs. 
Of course, because PAHs and grains are both modeled 
the same way, i.e., as spherical conductors, all of the 
equations below would still be valid were we to replace 
PAHs with grains. 



Standard metal abundance. — As stated above, we assume 
for this case that all ions are HCO + molecules and neglect 
M + . We may solve for x e and x\ = Xh~co+ m * ne limits of 
low and high PAH abundance. In t he l imit o f low xpah, 
the rightmost terms in Equations (16) and (17) can be 
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€pah 
10" 5 10" 2 



ignored, yielding in steady state: 




-lag 



•EPAH 

Fig. 5. — Average charge state of PAHs as a function of PAH 
abundance (solid diamonds, left axis). Dashed lines show simu- 
lation results for fractional electron abundance x c (solid circles, 
right axis) and fractional ion abundance x\ (open squares, right 
axis). The shaded region marks observationally inferred PAH 
abundances, measured by number either relative to H2 (a:pAH> 
bottom axis) or relative to the PAH abundance in the diffuse ISM 
(depletion factor epAH = x 'pah/10~ 6 , top axis). Parameters for 
this run are a = 3 AU, = 10 -8 , S = 0.3 g cm -2 , and e gra i n = 
(kept at zero to isolate the effect of PAHs). The shift to (zTpah) = 
occurs when there are so many PAHs that they begin to adsorb 
most of the free charge. At this point Xi and x e diverge; see also 
Figure [7] 
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Fig. 6. — Same as Figure[5]but for grains. Parameters 
are a = 3 AU, = 10 — , S = 0.3 g cm -2 , and epAh — . «• ,.. 
at zero to isolate the effect of grains). Our grains all have radii 
of 1 /xm, which is so large that their corresponding abundance as 
inferred from observation (shaded region) is too low to significantly 
affect the amount of free charge. 



X c = XHCO+ = V C/"H 2 aHCO+,, 

L x 



(18) 



10 



-in 



10 29 ergs" 



1/2 



) (iio) < 19 > 



v 0.3 gem 2 / 
for standard metals and low PAHs. 



In going from Equations (18 1 to (19) we account for the 



distance dependence of temperature but assume material 
is opt ical ly thin to X-rays. The square-root law of Equa- 



tion Jl8l) is often used by other workers (e.g., Gammie 



1996 jGlassgold et al.|1997[ ). It is plotted as a horizontal 



dashed line in Figure |7J and should be compared with 
the curves for x c and from our code, plotted as solid 
lines. In the limit of low xpahj the electron and ion 
abundances computed from the code are nearly constant 
with xpah, as predicted by the analytic model. However, 
the results from the code sit above the line for Equation 
( 18 ) by about an order of magnitude. The factor of 10 
offset arises because Equation ( 18 I ignores ionized met- 
als, which recombine with electrons much more slowly 
than does HCO + and which remain abundant compared 
to HCO + in our standard model. The offset also implies 
that reducing the total metal abundance below that of 
our standard model (xu = 10~ 8 ) can only decrease x e 
and x\ by at most a factor of ~10. In this sense our 



uncertainty in the metal abundance (Section 2.4.1) has 
only a limited impact on the io nizatio n fraction, assum- 
ing xm < 10 -8 . See also Section 3.3.1 where we consider 
the case xu = 0. 

In the limit of high PAH abundance, electron recombi- 
nation on PAHs dominates electron recombination with 
HCO + . Low equilibrium abundances of free electron s im - 
ply the average charge on PAHs (Z) (Section |3.1[ ). 
Then the steady-state solutions to Equations ( 16 ) and 
(17) are, respectively, 



c 



x wo+ 



^PAH^H 2 ( Q! PAH,c)(Z)=0 

: c 

^PAH^H 2 (apAH,HCO+)(Z)=0 

for high PAHs. 



(20a) 
(20b) 



Note that in this limit of high PAH abundance, x c < 
xnco+ because «pah,c > o;pah,hco+; electrons move 
faster than ions. The remaining negative charge required 
to maintain charge neutrality is carried by PAHs. That 
x c =/= £hco+ is r elevant for the computation of Am and 
Re (Section |3.3| because Am depends on Xjjco+ (ions 
carry the bulk of the momentum in a plasma) while Re 
depends on x c (elec trons are t he m ost mobile charge car- 
riers) . Expressions (|20a| and (|20b[ ) are plotted as diago- 
nal dashed lines in Figure [7j they compare well with the 
full numerical results, shown as solid lines, in the limit 
of high xpah- 

The critical PAH abundance dividing these limits is 
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estimated by equating Equations (20a) to (18 1 



charges is 



x pah — 



C a HCO+,c 



lH 2 (t»PAH,e) 2 



(21) 



-5 x 10" 



(Z)=0 



10 29 ergs" 1 

-1/2 



1/2 



, 0.3 gem" 2 / V3AU 
for standard metals 



The critical value xpAH m & r ks the abundance at which 
PAHs start to reduce significantly the number of free 
charges, i.e., the abundance at which electron recombi- 
nation on PAHs becomes competitive with electron re- 
combination with molecular ions. It is plotted in Figure 
[7] as a vertical line and does reasonably well at delin- 
eating the regime where the ionization fraction does not 
depend on PAHs from the regime where it does. Note 
that possible PAH a bundances as inferred from observa- 
tions (Section 2.4.3) happen to straddle Xp AH . 




Fig. 7. — Ionization fraction as a function of PAH abundance for 
o;m = 10~ 8 (standard metal abundance), a = 3 AU, and £ = 0.3 g 
cm -2 . Dashed lines: asy mptotic values for x c and i HCO + of the 
simplified model of Section |3.2.1| Solid lines: simulation results for 
fractional electron abundance x c (solid circles) and fractional ion 
abundance x^,+ +iHp n+ ( open squares). The dotted vertical line 
marks xi „ (Equation 12 1 1 , which roughly divides the regime of 
"low PAH abundance" where electron and ion densities are equal 
and insensitive to PAH abundance, from the regime of "high PAH 
abundance" where the ion density exceeds that of electrons and 
both decrease approximately as I/^pah- The behavior at high 
PAH abundance is independent of the metal abundance; compare 
with Figure [8] The shaded region marks observationally inferred 
PAH abundances and happens to span the transition from low to 
high PAH regimes. Simulation data use e gra in = 10" 3 for epAH < 



If] -5 - f = 
e PAH > 10- 2 . 



10 for epAH 



io- 



and 



1(T 



for 



Metal-rich case. — Analogous results are obtained for the 
metal-rich case as shown in Figure [8j with the only dif- 
ference that M + replaces HCO + as the dominant ion. In 
the limit of low PAH abundance, the abundance of free 



x c = x M + = y / C/«H 2 a M +, l 



^3 x 10" 



10 29 ergs" 1 

-1/2 



1/2 



0.3 gem' 



3 AU 



-0.5 



(22) 



(23) 



for high metals and low PAHs. 



In going from Equations (22) to (23) we account for the 



distance dependence of temperature but assume material 
is optically thin to stellar X-rays. Just as assuming all 
ions took the form of HCO + in the standard model gave 
a lower limit (Equation |18[ ) for the ionization fraction, 
assuming that all ions take the form of metals gives an 
upper limit (Equation 22 ) because fast recombination of 
electrons with HCO + is neglected. 

The analogous asymptotic solutions in the high PAH 
limit are practically unchanged from Equations ( 20a ) and 
( 20b ) because the charging rates of PAHs by HCvJ^ and 
lvl + are similar; the mass of the HCO + molecule and that 
of a metal ion like Mg + are similar. The critical PAH 
abundance at which PAHs begin to reduce the number 
of free charges is 



C"M-t 



PAH ■ 



"-H 2 ( a PAH,e)^) =0 



(24) 
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for high metals 
and is confirmed by the code. 
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Fig. 8. — Same as Figure [7] but for the metal-rich case (% = 
10 -6 ). The curves for x e and x\ at high PAH a bun dance (oipah > 
ip AE , where &i>, H is now given by Equation [24| are essentially 
the same as in Figure]?] when PAHs dominate charge balance, the 
metal abundance ceases to matter. 
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3.3. Degree of Magnetic Coupling: Re and Am 

We measure the extent of the MRI-active column by 
means of the magnetic Reynolds number Re (Equation!!]) 
and the ion-neutr al c ollisional frequency Am (Equation 
[2]). Figures [9] and 10 show both dimensionless numbers 
as a function of the surface density £ = N [i penetrated 
by X-rays at a = 3 and 30 AU, respectively, over the 
range of observationally inferred PAH abundances. We 
overplot for comparison the solution obtained when we 
omit PAHs completely. 

In both Figures [9] and [lOj the middle panels display 
results for our standard model parameters: im = 10~ 8 
per H2, Lx — 10 29 erg s" 1 , and £cr = 0. In each of the 
panels on the left and on the right, we vary one of these 
parameters. We describe here results for our standard 
model and compare with other test cases in Section |3.3.1 



The middle top panels of Figures [9] and [T0| show that 



- 2 ), 
-10 4 



Re were the only discriminant, MRI-active surface layers 
could well exist, even with PAHs present. At surface den- 
sities £ ~ 0.3 g cm' 2 (column densities N ~ 10 23 cm" 
Re lies c omf ortably above the critical values of 10 2 - 
(Section 1.2 1 required for plasma to couple to the mag- 
netic field, tor a wide range of possible PAH abundances. 
If the critical Re ~ 10 2 (as assumed by BG), and if PAHs 
are at their lowest possible abundance as inferred from 
observation (cpah = 10~ 5 relative to the ISM; inverted 
triangles) , then the MRI-active layer could extend as far 
as £ ~ 20 g cm" 2 — if ohmic dissipation were the only 
limiting factor for the MRI. 

But ohmic dissipation is not the only factor. The same 
margin of safety enjoyed by Re does not at all apply to 
the ambipolar diffusion number Am, for any surface den- 
sity. Even in the unrealistic case that there are no PAHs, 
Am stays < 10 in the middle bottom panels of Figures [9] 
and 10 By comparison, values of A m exceeding 10 2 are 
reported by Hawley & Stone ( 1998|) as necessary for the 



MRI to excite turbulence in predominantly neutral gas 
When PAHs are present, Am barely exceeds 1, and then 
only for the low end of possible PAH abundances. Com- 
pared with ohmic dissipation, ambipolar diffusion seems 
the much greater concern for the viability of the MRI in 
disk surface layers. 

Values of Am{T,) and i?e(£) vary only slightly as the 
stellocentric distance increases from a = 3 AU (Figure 
9l middle) to 30 AU (Figure nOl middle); the former de- 







creases while the latter increases, each typically by fac- 
tors of a few. This behavior is readily understood. First 
recognize that xi oc a~ o e approximately; this is an aver- 
age scaling between the low PAH li mit , which implies 
Xi oc a -05 according to Equation (18), and the high 
PAH limit, w hich implies Xi oc a~ 5 / 7 rs a~ ' 7 according 
to Equation ( 20b ) . Combining this result with njj 2 oc 
h- 1 oc O/T 1 / 2 ^- 



9 / 7 , we find that Am = x\nYi 2 /Sl oc 



-0.4 



Similarly, Re = c s h/D oc x^^/fl 



oc a 



0.7 



In computing Am, we have omitted the contribu- 
tion from collisions between neutral H 2 and negatively 
charged PAHs. The latter are as well coupled to mag- 
netic fields as molecular ions are — see, e.g., the ion and 
grain Hall parameters calculated in Section 2.2 of BG. 
Thus, collisions between H 2 and charged PAHs should 
increase Am. However, we find that in practice the gain 
is negligible. We estimate that the collisional rate co- 



efficient j3- m that enters into Am is about the same for 
charged PAHs as for ions; in both cases a collision with 
an H 2 molecule is mediated by the induced dipole in H2, 
and the relative velocity is dominated by the thermal 
speed of H 2 . At £ > 10 g cm~ 2 , charged PAHs are 
about as abundant as ions and thus raise Am by a factor 
of 2 — but at these £'s, Am is already too low for the 
MRI to be viable. At £ < 10 g cm -2 — i.e., at those 
columns where Am peaks — charged PAHs are much less 
abundant than ions and thus hardly affect Am. 

3.3.1. Higher Lx and Tx, Higher and Lower xm, and 
Cosmic-ray Ionization 

In each of the leftmost and rightmost panels of Figures 
[9] and |10| we vary one model parameter away from its 
standara value. We begin with the case of higher Lx- 
Increasing Lx certainly raises Re and Am, but as the left 
panels of Figure M show, even a fairly high Lx = 10 31 erg 
s _1 only causes Am to just exceed 10 at the lowest PAH 
abundance. In the limit of low PA H a bundance, x c — 

1/2 . . J L 

X{ oc L x , as predicted by Equation (18). In the limit of 

high PAH abundance, x 9 an d x\ s cale linearly with Lx, 
according to Equations ( 20a ) and ( 20b ) . Thus the space 
of possible values of Re and Am narrows with increasing 
Lx, as the lower envelope increases as Lx while the upper 

envelope increases as L-\[ 2 . 

These same scaling relations, with Lx replaced by the 
ionization rate £ at fixed distance, enable us to estimate 
the effects of a higher Tx, a case we did not explicitly 
compute using our numerical model. According to IG, 
raising kT x from 3 keV (our standard value) to 8 keV 
increases the ionization rate £ by factors of 2-4 at £ = 1- 
30 g cm -2 . Thus in the extreme case that Lx = 10 31 



erg 



s and kTx = 8 keV — parameter s appropriate only for 



Preibisch et al. 



a small minority of y oung stars (Telleschi et al. 2007 
12005) — we apply the low- PAH scaling 

m 



relation x\ oc c^ 1 / 2 to the bottom left panel of Figure 
to find that the largest possible value of Am is ^20, 
Dtained if PAHs are at their lowest plausible abundance. 
The rightmost panels of Figure M display the case of a 
higher me tal abu ndance xm = 10 per H 2 . As discussed 
in Section [2. 4. 1[ the higher metal abundance is not espe- 
cially realistic and is considered primarily as an exercise. 
Comparing the middle and right panels of Figure |9j we 
see that increasing the metal abundance by a factor of 
100 raises Re and Am at low PAH abundance by a fac- 
tor of ~10. At high PAH abundance, Re and Am also 
increase with increasing metal abundance, but the gain 
is less. This same behavior is reflected in the solid curves 
of Figures [7] and [HJ at low PAH abundance, increasing 
xm by a factor of 100 leads to a factor of ~10 increase in 
x c = Xi, but at high PAH abundance, the now divergent 
curves for x e and x\ are esse ntiall y ind ependent of metal 
abundance. Equations (20a) and (20b) from our analytic 



analysis reflect this insensitivity to metal abundance at 
high PAH abundance. 

At E > 1 g cm~ 2 , gas temperatures may be so low 
that all of the metals condense onto grains. The case 
xm = is shown in the leftmost panels of Figure [T0| 
Here Am < 0.1 for E > 1 g cm" 2 , and it seems safe to 
conclude that X-ray driven MRI is unviable under these 
conditions. 

Additional ionization by "sideways cosmic-rays" at a — 
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30 AU is considered in the rightmost panels of Figure 



10 These cosmic-rays, which we have imagined enter 
the disk edge-on from the outside, dominate stellar X- 
rays at large E. At S ~ 10 g cm -2 — comparable to the 
full surface density of the disk at a = 30 AU — sideways 
cosmic-rays raise the maximum value of Am to ~2. We 
have verified that the gains in Am afforded by cosmic- 
rays are consistent with our scalings of x e and x\ with £ 
as derived above. 

3.3.2. Chemical Equilibration Timescales vs. Dynamical 
Timescales 

In assessing whether disk surface layers are MRI-active, 
we have relied on the critical value Am* ~ 10 2 reported 
by HS. As a simplifying assumption, HS held fixed the 
global (box-integrated) ion abundance in each of their 
simulations. A fixed ion abundance would apply if the 
chemical equilibration timescale t cq exceeds the dynam- 
ical timescale <dyn = ovcr which HS's simulations 
ran. A fixed ion abundance would also apply if ion re- 
combination occurs pre dominantly on conde nsates, re- 
gardless of icqAdyn (e.g., Mac Low et al. 19 95). This last 
statement follows from our Equation ( 20b ) , which shows 
Xinu 2 does not depend on nn 2 in the high condensate 
limit. 

The high condensate limit applies for PAH abundances 
near the high end of those inferred from observation (Fig- 
ure [7]). For this high PAH case we expect the assumption 
of constant ion abundance, and by extension the results 
of HS, to hold. For high PAH abundance and our stan- 
dard X-ray luminosity, Am < 1 for all E and a (Figure 
[9]) , and our conclusion that X-ray-driven MRI shuts down 
everywhere seems safe. 

For PAH abundances at the low end of those inferred 



from observation, Equation (20b) for the high conden- 



sate limit does not apply. Moreover, as shown in Figure 
^cqAdyn < 1 — but only by a factor of 10 at most 



11 



for the low PAH case and for E < 10 g cm -2 . Because 
ieq/idyn ^ 0.1 under these conditions, the assumption 
of constant ion abundance, although not strictly valid, 
might still be good enough that the results of HS hold to 
order unity. But even if they do not, we would argue that 
the sign of any correction for a dynamically variable ion 
abundance would only hurt the prospect s for MRI turbu- 
lence. In the si mul ations of HS — see also Brandenburg fc| 
Zweibel| ( [19941 ) and |Mac Low et al.| ( p95| — ions became 
concentrated in thin filaments within magnetic nulls. Re- 
combination rates inside the dense filaments were higher 
than those outside. Were such simulations to account 
for ion recombination, lower ion densities within the fila- 
ments would result, and neutrals would be even less cou- 
pled to ions. Consequently, Am* would be even higher 
than the reported value of ~10 2 . 



3.4. Comparison with Previous Work: Ionization 
Fractions 

The ionization chemistry in disks remains inherently 
uncertain, with rate coefficients for many reactions in 
the UMIST database determined t o no better than fac- 
tors of ^3 (Vasyunin et al. 2008). A measure of the 



uncertainty in the ionization fraction in disks is given by 



the differences between the simple and complex networks 
computed by BG, which amount to factors of 2-10 for the 
electron fraction x c . We should reproduce their results 
by at least this margin, as a validation of our code. In 
the following we directly compare our results to those 
of BG and TCS, adjusting the input parameters of our 
code to match theirs. Once we match these input param- 
eters, any difference in our codes' outputs should result 
primarily from our different chemical networks (ours is 
the simplest of the three), and not from differences in 
radiative transfer, as all our codes rely on t he ionization 
rates calculated by Igea & Glassgold (19991. 

We start with BG by computing x e as a function 
of density njj 2 at a fixed ionization rate £ = 10 -17 
s" 1 , following their Figure 3. We reset T = 280 K, 
im = 2.5 x 10~ 8 per H2, and the electron-grain sticking 
coefficient S e — 0.03 to match their standard parame- 
ters. To compare to their "grain-free" case, we run our 
code without any PAHs or grains. To compare to their 
standard monodispersion of grains, we run our code with 
a single population of grains having s = 0.1 pm, internal 
density p s — 3 g c m - 3 , and a mass fraction of 1% rela- 
tive to gas. Figure 12 shows the comparison. Our results 
for the condensateTree case track those of BG, but are 
higher by factors of 3-10 depending on whether the com- 
parison is made with their simple or complex network. 
For the case with grains, the agreement with the simple 
model is excellent and that with the complex model is 
good to a factor of 2. 

In Figure 13 we make a similar comparison with TCS, 
computing x c as a function of A at a distance of a = 5 
AU from an X-ray source of Lx = 2 x 10 30 erg s _1 and 
kTx = 5 kcV, for T = 125 K and a metal abundance 
of xm = 6.8 x 10 -7 per H2. We consider the two cases 
of their Figure 1, one without any grains or PAHs, and 
another with a single population of grains having s — 
1/xm, p s = 5 g cm~ 3 , and a mass fraction of 1%. For 
both cases our computed electron abundances are higher, 
but only by factors of 2 or less. 

These comparisons with BG and TCS give us confi- 
dence that we have computed ionization fractions about 
as well as they did. Where our ionization fractions differ, 
ours are often higher. Our higher values will only bolster 
the conclusion we make in Section 2] that thicknesses of 
X-ray-ionized MRI-active surface layers have been over- 
estimated by them and others. 

4. SUMMARY AND DISCUSSION 

In Section [T] we presented the evidence that holes and 
gaps of transitional disks are cleared by companions to 
their host stars. Residing within the hole, these compan- 
ions could either be stars — already observed in about half 
of all transitional systems — or multi-planet systems. A 
single Jupiter-mass planet on a circular orbit carves out 
too narrow a gap to explain the large cavities inferred 
from observations. But multiple planets can shuttle gas 
quickly from one planet to the next, all the way down to 
the central star. Surface densities fall in inverse propor- 
tion to radial infall speeds, and radial infall speeds can 
approach freefall speeds for sufficiently many and mas- 
sive planets. In this way, multi-planet systems might 
help to clear holes and simultaneously sustain stellar ac- 
cretion rates that approach those in disks without holes. 
The more eccentric the planets' orbits, the fewer of them 
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Fig. 9. — Magnetic Reynolds number Re and ambipolar diffusion number Am as a function of surface density S at a = 3 AU. The 
middle panels show results for our standard model (xm = 10 — 8 , Lx = 10 29 erg s _1 , Ccb. = 0). The side panels have the same parameters 
as our standard model, except for a lOOx more luminous X-ray source (left panels), and a lOOx greater metal abundance (right panels). 
Oppositely point ing triangles bracket values for Am and Re corresponding to possible PAH abundances. These abundances were inferred 
in Section |2.4.3| from observations. The dashed curve refers to the case with no PAHs and is shown for comparison only. Polycyclic 
aromatic hydrocarbons reduce ionization fractions and thus the degree of magnetic coupling by an order of magnitude or more. The dotted 
lines ma rk t he critical values Am* and Re* above which coupling between magnetic fields and neutral gas is sufficient to drive the MRI 
(Section|P). The curves for Am first rise as £ increases — a consequence of the increasing number density — and then fall as the ion fraction 
decreases, never reaching Am*. Ambipolar diffusion threatens the MRI more than ohmic dissip ation does. If the critical Am required for 
good collisional coupling between ions and neutrals is Am* = 10 2 , as evidenced in simulations by |Haw!ey fc Stone| ( [l998| , then even a disk 
without PAHs cannot sustain X-ray-driven MRI at any S. 



PAHs, MRI, and Planets 



19 



may be required to explain a given hole size. Accretion 
in the presence of multi-planet systems has not received 
much attention and seems an interesting area for future 
simulation (e.g., Zhu et al. 2010, submitted). 

Stellar or planetary companions regulate accretion ve- 
locities v but do not give rise to mass accretion rates M 
in the first place. A planet orbiting just inside the cir- 
cumference of a disk hole exerts torques to repel gas in 
the hole's rim away from the star. Thus, the shepherd- 
ing planet may reduce M — and indeed accretion rates 
in transitional syst ems tend to be sma ller than those in 
conventional disks (Najita et al. 2007) — but the planet 
does not initiate disk accretion. A separate mechanism 
must act to pull or diffuse gas inward from the hole rim to 
supply the stellar accretion rates that are observed. That 
mechanism may be turbulence driven by the magnetoro- 
tational instability (MRI), activated by stellar radiation 
ionizing rim gas. Whether the MRI can operate depends 
on how well ionized the gas is. The greater the free elec- 
tron fraction, the greater the magnetic Reynolds number 
Re, and the less ohmic dissipation dampens the MRI. 
The greater the atomic and molecular ion densities, the 
greater the collisional rate Am between neutral particles 
and ions, and the less ambipolar diffusion weakens the 
MRI. 

A principal threat to the MRI is posed by dust grains, 
which adsorb electrons and ions. The smallest grains 
may present the biggest danger, because in many parti- 
cle size distributions the smallest grains have the great- 
est surface area for attachment. The smallest grains 
that can also be detected observationally are polycyclic 
aromatic hydrocarbons (PAHs), each several angstroms 
across and containing of order a hundred carbon atoms. 
These macromolecules may reside in the very disk sur- 
face layers that promise to be MRI-active. Excited by 
soft ultraviolet radiation, PAHs fluoresce in a distinctive 
set of infrared emission lines detectable from Spitzer and 
from the ground. The hydrocarbon molecules are proba- 
bly generated locally, photo-sputtered off larger particles 
exposed to hard UV and X-ray radiation from host stars. 

To assess the impact of PAHs on the MRI, we need 
to know PAH abundances relative to gas. These can 
be inferred from observed PAH emission lines. Unfortu- 
nately such inferences are model dependent; they depend 
on knowing the local grain opacity, because the soft UV 
radiation that excites PAHs is also absorbed by grains. 
In other words, observed PAH line intensities depend on 
PAH-to-dust ratios. It follows that the quantity of inter- 
est to us — the PAH-to-gas ratio — depends on knowing 
the dust-to-gas ratio. The latter can vary widely with 
the degree to which grains settle toward disk midplancs. 
The more grains have settled, the lower are local dust-to- 
gas ratios, and the lower the PAH-to-gas abundance that 
is needed to explain a given set of PAH emission spectra. 
By compiling a few lines of model-dependent evidence 
from the literature, we estimated that disk PAHs have 
abundances anywhere from 10 _11 -10~ 8 per H2, with 
lower values corresponding to a greater degree of dust 
settling. 

Such PAH abundances, although depleted relative to 
the ISM by 10 _5 -10 -2 , are still large enough to signif- 
icantly weaken the MRI in disk surface layers. In fact, 
they might even shut off X-ray-driven MRI altogether, 



everywhere. For stellar X-ray luminosities of L> 



10 29 - 

10 31 erg s _i and X-ray stopping columns of £ ~ 1-10 g 
cm~ 2 , PAHs reduce electron and ion densities — which are 
not equal when PAHs are present — by factors of ~10 or 
more. At these surface densities, the collisional coupling 
frequency Am ss 10 _3 -10, depending on PAH abundance 
and X-ray luminosity. These values fall short, by 1-5 
orders of magnitude, of the critical value Am* ~ 10 2 re- 
quired for good coupling b etween ions and neutra ls, as 
measured in simulations by Hawley & Stone (1998) HS). 
The potentially catastrophic effect that small grams can 
have on the MRI was highlighted by Bai & Goodman 
(2009, BG). Our study grounds their concern in real-life 
observations. 

Other studies reported X-ray-driven MRI-active sur- 
face layers to be alive and well (e.g., Chiang & Murray- 
Clay 2007, CMC; Turner, Carballido, & Sano 2010, TCS; 
and BG, in many of whose models the active layer ex- 
tended to ~1 g cm~ 2 , even with grains present). We 
should understand why our conclusions differ from theirs. 
In part, the difference arises because previous studies ne- 
glected PAHs. A further difference with TCS is that they 
assumed a metal abundance of xm = 6.8 x 10 -7 per H 2 , 
nearly 2 orders of magnitude higher than our standard 
model value, and one that we find difficult to justify. Still 
another difference, as significant as any of the ones just 
mentioned, is the criterio n used for whether a mbipolar 
diffusion defeats the MRI. Turner et all (|2010|) assumed 



Am* — 1 (see their Equation 9). Bai fc Goodman ([2009) 
did not present results for Am. Using their data, we com- 
puted the Am values characterizing their claimed active 
layers. At a > 1 AU, the Am values of BG's grain-free 
active layer are at most on the order of unity. For BG's 
standard models containing grains, Am s» 0.0004-0.4, 
with the lowest value corresponding to a = 50 AU and 
a population of grains having two sizes, and the highest 
value corresponding to a = 1 AU and a single-sized grain 
population (we computed both li mits using results from 
their complex chemical network). Hawley & Stone ( 1998 1 
showed that when Am < 0.01, ions and neutrals were ef- 
fectively decoupled. Even when Am ~ 1, HS showed 
that the MRI saturation amplitude scaled with the ion 
and not the neutral density, with the neutrals acting to 
damp out MRI turbulence in the ions. If the MRI drives 
turbulence only in the ions of protoplanetary disks, it 
might as well not operate at all, given how overwhelm- 
ingly neutral such disks are. 

4.1. Future Directions 

We have shown in this paper that the MRI cannot 
drive surface layer accretion under typical circumstances 
in protoplanetary disks, either transitional or conven- 
tional, if the critical Am* ~ 10 2 and if stellar X-rays and 
Galactic cosmic-rays are the dominant source of ioniza- 
tion. These two "if" s are subject to further investigation. 
We discuss each in turn. 

We are not aware of more modern estimates of Am* 
apart from that given by HS. As these authors cautioned, 
numerical resolution is a greater concern for two-fluid 
simulations than for single-fluid ones, and HS did not 
demonstrate convergence of their results with resolution. 
In addition, the value of Am* was not as precisely deter- 
mined by HS for toroidal field geometries as for vertical 
ones — although Am* ~ 10 2 did seem to apply equally 
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Fig. 10. — Same as Figure [9] but at a stellocentric distance of a = 30 AU. The middle panels show results for our standard model 
(iM = 10 — 8 , Lx = 10 29 erg s _1 , Ccr = 0). The side panels have the same parameters as our standard model, except that gas-phase 
metals are omitted in the left panels (m = 0), and sideways cosmic-rays are added in the right panels (Ccr = 1/4 X 10 — 17 s — 1 ). 
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Fig. 11. — Ratio of the chemical equilibration timesca le t ca , com- 
puted according to the procedure described in Section |2,6| to the 
dynamical time t^ yn = S7 — 1 , for a = 3 AU and = 10 . Values 



of t e q/tdyn at 



30 AU are typically lower than those shown 



here by factors of 3 or less. Simulations by HS assumed a con- 
stant (volume-integrated) abundance of ions, a condition satisfied 
if teq/idyn > 1- The ion abundance is also constant if ion re- 
combination occurs primarily on condensates (see Equation |20b| , 
a situation that obtains for PAH abundances near the high end 
of those inferred from observations. For the lowest possible PAH 
abundances, 0.1 < teq/*dyn < 1 at £ = 0.1—10 g cm -2 . In this 
low PAH case, the results of HS might still be expected to apply 
to order unity. Even if they do not, we argue in the main text that 
when icqAdyn < 1 the effects of ambipolar diffusion should be even 
stronger than reported by HS. 
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Test comparison with BG 
(fixed C = 10- 17 s- 1 H^ 1 ) 
Our code 
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Fig. 12. — Test comparison with BG: electron abundance as a 
function of gas density at a fixed ionization rate of £ = 10~ 17 
s" 1 H" 1 . The upper set of lines are for a condensate- free system, 
and the lower set are for a monodispersion of grains. Results from 
BG were drawn from their Figure 3. To generate our results, the 
parameters of our code were reset to those of BG: temperature T = 
280 K, metal abundance ijj = 2.5 X 10 -8 per H2, electron-grain 
sticking coefficient S c = 0.03, grain radius s = 0.1 fim, internal 
grain density p s = 3 g cm -3 , and a mass fraction in grains relative 
to gas of 1%. 



well to the cases of uniform vertical field and zero net 



Fig. 13. — Test comparison with TCS: electron abundance as 
a function of column depth penetrated by X-rays. The upper 
pair of lines are for a condensate-free system, and the lower pair 
are for a monodispersion of grains. Ionization rates are from 
IG. Results from TCS are taken from their Figure 1. To gen- 
erate our results, we reset the parameters of our code to match 
those of TCS: stellocentric distance a = 5 AU, X-ray luminosity 
Lx = 2 X 10 30 erg s , temperature T = 125 K, metal (magne- 
sium) abundance xm = 6.8 X 10 — 7 per H2, electron-grain sticking 
coefficient S B = 0.03, grain radius s = 1 /im, internal grain density 
p s = 5 g cm -3 , and a mass fraction in grains relative to gas of 1%. 

vertical field. Perhaps higher resolution simulations will 
reveal that Am* < 10 2 — although accounting for ion re- 
combination i n the se simulations should only increase 
Am* (Section [3X21. 



The second possibility is that our model has neglected 
a significant source of ionization. Stellar radiation just 
longward of the Lyman limit — so-called far ultraviolet 
(FUV) radiation at photon energies between ~6 and 13.6 
eV — can ionize trace specie s such as C, S, Mg, Si, and Al 
(e.g., Cruddace et al.(l974 |. Of these elements, C and S 
may be the most important, as they are the most abun- 
dant and possibly least depleted onto grains. For exam- 
ple, if the full solar abundance of C were singly ionized, 
the ion fraction x\ would be a few x 10 -4 , or 5 orders of 
magnitude higher than the largest values of x\ reported 
in this paper! At disk midplanes which are shielded from 
photodissociating radiation, an order-unity fraction of 
the full solar abundance of C is expected to take the 
form of CO (xco = 1CT 4 ; Aika wa et al.||l996). A s com- 
puted i n chemical models by Gorti &; Mollenbach (2004; 
see a lso Ticlens & Hollenbach 19851 and IKaufman et al.l 



19991, CO near disk surfaces photodissociates nearly en- 
tirely by FUV radiation into a layer of neutral C. At the 
highest altitudes, nearly all of this carbon is photoionized 
by FUV radiation. The column density of C + depends on 
how many small grains having sizes < 0.1/im are present, 
as grains compete to absorb the same FUV photons that 
photodissociate CO and photoionize C. 

We may estimate maximum FUV-ionized column den- 
sities by neglecting such dust extinction, and by neglect- 
ing shielding of FUV radiation by molecular hydrogen. 
Consider a trace species T whose total number density 
regardless of ionization state is /t^h 2 • Take all of T to 
be singly ionized within a Stromgren slab at the disk sur- 
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face: nx+ = n e = /t^h 2 - P er un it surface area of slab, 
the rate of photoionizations balances the rate of radiative 
recombinations : 



In Equation (27), we have normalized /t to its highest 



J FUV 



E-puyATra 2 



(25) 



where the FUV lumi nosity capable o f ionizing T is 
Lpuv ~ 10 30 erg s" 1 ( | Gorti et al.|2009[ ), the photon en- 
ergy £fuv ~ 10 eV, the rate coefficient «T+,e ~ 4x 10~ 12 
cm 3 s _1 at an FUV- heated gas temperature of 300 K, 



and the slab thickness h 
column 



0.1a. Solve for the hydrogen 



A i 



FUV : 



-m 2 h 



■2 x W 



(26) 



( L 



FUV 



3AU 



VlO 30 

1/2 



ergs 
7r~ 



1/2 



or equivalcntly 



^fuv = Npuvll 



(27) 



-0.07 



( L 



FUV 



1/2 



\10 30 ergs" 



3AU\ 1/2 /10" 4 



h 



gem 



plausible value, appropriate for C. An MRI-active surface 
density Spuv ~ 0.07 g cm" 2 is modest, and would drive 
mass accretion rates only barely observable. Lowering 
/t would increase Spuv- But accounting for extinction 
of FUV radiation by dust and molecular hydrogen would 
decrease Spuv- We are currently undertaking a more 
careful study of FUV ionization to quantify these effects. 
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Note added in proof. — In a private communication, Xue-Ning Bai reports that in new, unstratificd MRI simulations, 
the Shakura-Sunyaev transport parameter a is at most ^10~ 3 when Am ~ 1. Combining his result with the results of 
our paper, we estimate that the mass accretion rate in the surface layer of a conventional disk at 3 AU is ^1O _11 M0 
yr _1 . At the rim of a transitional disk, we would predict an accretion rate that is lower by a factor of ^h/a, or 
^1CP 12 M Q yr _1 at 3 AU. These theoretical accretion rates are too low to explain the observed accretion rates of most 
disks. The situation is similar at 30 AU, where the surface layer accretion rate in conventional (transitional) disks can 
only be as high as ^10~ 9 (10~ 10 )M Q yr -1 if Galactic cosmic-rays can penetrate the disk edge-on (see our Figure 10, 
right). 



